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We apply spin-squeezing techniques to identify and quantify highly multi-partite photonic entan- 
glement in polarization-squeezed light. We consider a practical single-mode scenario, and find that 
Wineland-criterion polarization squeezing implies entanglement of a macroscopic fraction of the to- 
tal photons. A Glauber-theory computation of the observable A-photon density matrix, with A up 
to 100, finds that A-partite entanglement is observable despite losses and without post-selection. 
We estimate that existing detectors could observe ~ 1000-partite entanglement from a few dB of 
polarization squeezing. 



Introduction - Multi-particle entangled states are of cen- 
tral importance in quantum computing and quantum 
simulation, where entanglement gives an advantage that 
grows rapidly with the size of the entangled system. 
Multi-particle entanglement has been demonstrated by 
controlled interactions of individual particles in a vari- 
ety of systems, including trapped ions [1^, liquid-state 
NMR [2], solid state devices and photons [IHS]. The 
current record is 14 particles [7| (ions). 

A complementary approach works with large en- 
sembles and uses squeezing, either dynamical or 
measurement-based, to generate a non-classical state of 
the ensemble. S0rensen and coworkers |8] showed that 
for spin-1/2 atoms, spin squeezing implies entanglement 
within the ensemble. Similar results have been found 
for higher-spin atoms PJ, for multi-partite entanglement 
[TUHIS s-'^d fo'^ other kinds of squeezing [O]. Squeez- 
ing can generate highly multi-partite entanglement, the 
"extreme spin squeezing" of [10] , because it acts simulta- 
neously on a large number of particles. Squeezing of spin 
ensembles has been reported |14H19| . implying pairwise 
entanglement of macroscopic numbers of particles and at 
least 80-partite entanglement [14], although microscopic 
correlations have not been directly observed. 

Here we consider the analogous situation for photons. 
As photons have two possible polarizations, the natural 
analogy is to squeezing of spin-1/2 ensembles. Unlike 
atoms, no superselection rule prevents photonic systems 
from existing in a superposition of different particle num- 
bers. Indeed, these superpositions are the defining fea- 
ture of optical coherence [201 HI]- Because optical co- 
herence at all orders is abundantly present in laser light 
[20j . it is arguably a classical, not a quantum, resource. 
Optical coherence is nonetheless very useful for gener- 
ating photonic entanglement, as we show. We describe 
a simple scenario that requires as input only one weak 
squeezed- vacuum state and one coherent state, and which 
generates A^-partite entanglement for arbitrarily large N. 
We estimate that ~ 1000-partite entanglement should be 
directly observable with existing detectors. 

Scenario - We consider a single-mode pulse of 
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FIG. 1. Scenario for entanglement generation: a horizontally- 
polarized coherent state Ian) is combined with a vertically- 
polarized squeezed state S{r)pthS (r). The resulting 
polarization-squeezed state is symmetrically dispersed by a 
lens to an array of polarization analyzers (PA), each consisting 
of wave-plates, a polarizing beamsplitter, and single-photon 
detectors. Under conditions typical of squeezing experiments, 
nearly all of the photons originate in the coherent beam and 
become polarization entangled via coherence with the much 
weaker squeezed beam. 



polarization-squeezed light 



described by the state 
(1) 



P= Ph® Pv 



consisting of a coherent state pu = |a//)(ai/| in the 
horizontal [H) polarization and a squeezed thermal 
state Pv = S{r)pt\^S'^{r) in the vertical (V). Here 
S\r)avS{r) = aycoshr- aySinhr defines the squeeze 

operator S and pth = Zn'^"h/(1 "*" ^^th)"^^ |^) (j^l is the 
thermal state. It is convenient to parametrize this state 
by He = |ap. Us = sinh \r\ and nth- States with this de- 
scription can be produced by combining squeezed and 
coherent pulses on a polarizing beam-splitter [24] or by 
polarization self-rotation [25] of a coherent pulse. Re- 
garding py, for simplicity we take r real and negative 
and define T^ = 1 + 2nth, A = ^ri^ + \/ns + 1, so that the 
quadratures x = {ay + ay)/2, p = i(ay - ay)/2 have vari- 
ances var(a;) = T^A/4, var(p) = T^/iiA^). Note that p is 
squeezed if Ug > n^^^Kl + 2nth)- 

Macroscopic quantum features - The quantum description 
of the macroscopic polarization is given by the Stokes 
operators [55] Sq = {a\jaH + CLyO-v)l'2, Sx = (a'^an - 
a[-ay)/2, Sz = -i{a^^av -a|^a/f)/2, iSy = [Sz^S^]- We 
find 2 (5o) = n^ + n^ + ^^th + Srignth, 2 {S^) ^n^-Us- nth " 
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FIG. 2. (color online) Entanglement depth as a fraction of 
total photon number fc/(2 (So)) in the thermodynamic limit 
(So) -> oo, shown versus squeezing n^ and thermal noise nth- 
Contour values shown on right. Grey region indicates non- 
entangled regime k = 1, which coincides with non-squeezed 
states. 



2?^s"-th- If a = \ol\ so that S^ is the lower-noise Stokes com- 
ponent, we find 4var(S'2) = nc{\ + 2ns-2^ns{ns + !))(! + 
2nth). The state is squeezed by the Wineland criterion 



71s > n 



th 



/(271th + l)i the same 



EH, var(5,)<i|(5,)|if 
condition as for pv quadrature squeezing. 
Microscopic quantum features - To find the depth of en- 
tanglement fc, i.e., the minimum size of entangled clus- 
ters in the state, we use the theory of S0rensen and 
M0lmer [TU], extended by Hyllus et al. [21] to indefinite 
numbers of particles. For given degree of polarization 
V = {Sx) / (So) and polarization noise V = var{Sz)/ {So), 
k/2 is the minimum spin capable of producing {V,V). 
As 71c ^ oo, V ^ 1(1 + 271s - 2 \/ns(7is + !))(! + 2nth) 
while V ^ I. If V < |, i.e., if the state is squeezed, 
then P ->■ 1 implies k ^ oo, indicating of arbitrarily large 
groups of entangled particles. As shown in Appendix 
I, in this limit k is proportional to the total number of 
photons 2{Sq). Polarization squeezing evidently implies 
a giant connected component in the language of graph 
theory PU and entanglement percolation |30J EI] . The 
connected component can include nearly all photons for 
modest rig, iith, see Figure [2] 

Observable entanglement - Reconstruction of density ma- 
trices from coincidence rates is the standard method for 
characterizing multi-photon states |32| . The coincidence 
rates are described by Glauber photo-detection theory 
[20] , while the reconstruction employs quantum state to- 
mography. The results indicate not the ideal p of Eq. ([l} , 
but rather an observable density matrix (ODM) TZ^ ' for 
each possible number of detected photons N. This prop- 
erly accounts for loss of information to the environment 
in the form of undetected photons; TZ^ ' is usually mixed 
even when p is a pure state. 

Concretely, we imagine passively splitting the beam 
to a large number M » TV of analyzers, as illustrated 



in Figure [T] and using the observed detections to in- 
fer the locations of the photons. The d'th detection 
mode is described by the annihilation operator a^^p = 

where p is the polarization and 



VW-^Qp+VT- 

(aux) 



M 



1 (aux) 
d,p 



a^ describes an auxiliary mode, required for unitarity 
of the splitting process and assumed to be in a vacuum 
state. Because the splitting is polarization-independent, 
the output retains the exchange symmetry of the initial 
state p, and thus only the polarizations of the detected 
photons, not their arrival locations, are relevant to the 
entanglement characterization. 

For a single pulse, the probability of coincidence detec- 
tion with polarizations p^^^ , . . . , p^^^ at detectors 1 . . . A^ 
is given by 

d(A') / t t \ 

°^ ("p(i) • ■ • "p(") V"> • • ■ «p(i) / ' (2) 

where the second line follows from the definition of 
ad^p and the fact that the vacuum contributes noth- 
ing to the normally-ordered expectation. We employ 
the computational basis Ittq) = \H, ...,11, H) , \-Ki) = 
\H,...,H,V),..., |7r2iv_i) = \V,...,V,V) to describe the 
(unnormalized) ODM 
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), (3) 



where tt^ indicates the polarization of the Tth pho- 
ton in the basis element JTr^). We normalize 'RS'^' = 
i?(^)/Tr[i?(^)] so that 



Pj^/ = Tr[^p^7^(^)] 



(4) 



(0 



resembles the Born rule. Here p is defined to give p 
T,j Pj'^j , so that lip =p/\p is a projector onto {p}. 

Nature of the entanglement - TiS ' is in general a mixture 
of symmetric states, e.g. VS^' is a mixture of \H, H) ± 
\V,V), \H,V) + \V,H), and either \H,H) or \V,V). We 
show TiS^' and TiS^' as examples. 



7^(2) 



I A2 ■ ■ D2\ 
B2 B2 

B2 ■ 

\ h.c. C2 I 



(5) 



where the basis is \H, H) , \H, V) , \V, H) , \V, V) and A2 

37l2(l + 271th)' 



71 



^ B2 = nc{ns + 71th + 27is7ith), C2 
71s (1 + 87ith + 1271th) + 2'^th' and D2 



+ 



jij) + 2n'^^, and D2 = nc\/nJn^Tl) (1 + 
271th). The symbol "•" indicates a zero entry, while "h.c." 
indicates the lower-triangular portion of the (hermitian) 
matrix. In the limit of zero noise 7ith -*■ 0, and small flux 
71s = 71^ ^ 0, 7^(2) describes a Bell state oc \HH) + \VV). 



Similarly, 



7^(3) 



/ A3 ■ ■ E3 ■ E3 E3 ■ \ 

B3 B3 ■ B3 • • F3 

B3 • B3 • • F3 
C3 ■ C3 C3 

B3 ■ ■ F3 

h.c. C3 C3 ■ 

C3 ■ 



(6) 



where A3 



nl,B3 = nKris + nth + 2?isnth),C3 
nc[3n2(l + 2nth) + ns(l + 8nth + 12^^^) + 2n^^],D3 
3(ns + nth + 2r7,snth)[2njij + ns(l + 27ith)(3 + lOnth) 



5(7is + 2nsnth) ],E3 



?v/^I^(lT^(l + 2nth),-p3 



3nc\/ns{l + 7is)(l + 2nth)(ns + nth + 2nsnth)- In the same 
limit, 7e(3) describes a state oc \HHH)+\HVV) + \VHV) + 
\VVH) with maximal 3-tangle [55] , 
Multi-partite entanglement - To explore the entanglement 
properties with larger photon numbers, we study the re- 
duced two-photon density matrix, obtained by tracing 
over the last N -2 photons 



7^('''^)=Tr3...^^[7^(^)]. 



(7) 



By permutation symmetry, an entangled TZ^^' ' implies 
that no partition of TZ^^' is separable, i.e., A^-partite 
entanglement. As shown in the Appendix, TZ^"^'^' can be 
computed in 0{N^) time using phase-space methods. 
Reduction to pure states - We can simplify calcula- 
tion of TZS^' by considering the effect of polarization- 
independent losses. If Oi^TT -*• ai^Tr\/v + a^°^' ^yi - 77 
where rj is the transmission and a]^^^ describe reservoir 



modes (assumed vacuum), then R^ ' 



loss 



V 



N 



i?(^) and 



7?.^^ is unchanged. At the same time, Uc -^ rju^, and 

var(g) -^ r]var{q) + {I - r])/4, where q e {x,p}, from 
which we find that any mixed squeezed state p, with 
Uc, Us, nth, can be produced from a pure state p' with 
n'^ = nc/77, n^ = (ng + 2nsnth + nth)/?7 and n[^^ = for 
77 = (ns + 2nsnth - n^j^)/(ns + 2nsnth + nth)- 
Likely outcomes are entangled - For pure states, i.e., 
nth = 0, we compute the observable concurrence C, i.e., 
the concurrence computed on 7?."' % shown in Figure 
[3| In all cases, Wineland-criterion squeezing ng > 
gives nonzero concurrence. For example, with n^ = 100, 
Us = 0.3, and N = 100, 



/ 0.9440 



7^(2,100) 



0.0293 \ 



0.0270 
0.0270 



0.0270 
0.0270 



\, 0.0293 



0.0021 / 



(8) 



The concurrence is C « 0.00468, which is 4.7% of that 
allowed by monogamy [33], i.e., Cmax = 1/V-/V - 1 « 0.10. 
7^(2, Af) (Joes not depend strongly on N, and the likely 
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FIG. 3. (color online) Scaled concurrence C/Cms 
of A*', Uc, and n^ for pure states, i.e., rith = 
I/VN - 1 is the largest concurrence allowed by entanglement 
monogamy [34], shown for N = 2 (red), N = 17 (orange). A'' = 
50 (green) and A*' = 100 (blue). Labels indicate ris values, and 
line strength indicates the likelihood of observing N photons 
from the coherent state with n^ average photons, calculated 
as Pjv(nc)/-Pjv(A''), where Piv(A) is the Poisson distribution. 



outcomes Uc - y/n^ < N < n^ + ^J^\: have similar C/Cma^. 
In this sense, post-selection is not necessary to prepare or 
observe the entanglement. For these likely states, C peaks 
for Us ~ 0.3, or « 4.5 dB of squeezing. We conjecture that 
3-tangle and higher-order entanglements compete with C 
as ng grows. 

Prospects for observation- Due to its structure, TZS'^' ' is 
negative under partial transpose if A e \TCy ^ \~T^2 3 ^ 
0, e.g. if 0.0293 > 0.0270 in 7^(2.loo) above. For nc = iV in 
the range 2 to 100 and ng chosen to maximize C, we find 
A = 0.23/ N to within 10%. Although reconstruction of 
general A'^-particle states requires resources exponential 
in A'^, here the method of preparation ensures permuta- 
tion invariance (PI), allowing use of PI reconstruction 
methods to avoid exponential costs [3S] . In particular, 
we can find R^"^'^^ by averaging all pairs of detected pho- 
tons, i.e., N{N-1) pairs per shot, giving single-shot sta- 
tistical errors of ^A ~ 1/A^. We then need 0{N^) shots 
to establish A > with unit signal-to-noise. In light of 
the A^*^ scaling of both production and characterization 
of entanglement, A'^ may be limited by the performance 
of available detectors. In this regard, electron-multiplied 
charge-coupled device (EMCCD) sensors are promising. 
They can have M = 2^° > 10^ pixels, 90% quantum ef- 
ficiency, negligible dark-count rates, single-photon sensi- 
tivity, and some degree of photon-number discrimination 
for < 5 photons at a given pixel [3S] . Evaluating the effect 
of imperfect detection on TZ is beyond the scope of this 
work, but we note that for A^^ < M, fewer than one pixel 
on average receives multiple photons, so that entangled 
states with A^ ~ 1000 should be detectable even without 
photon-number resolution. 



Conclusion - In analogy with "extreme spin squeezing," 
in which squeezing of collective spin variables implies 
highly multi-partite entanglement of the constituent par- 
ticles, we have quantified the photonic entanglement of 
a practical polarization-squeezed state. We find that 
Wineland-criterion squeezing implies entanglement of a 
macroscopic fraction of the total photons present. It 
also implies A^-partite entanglement in the observable A^- 
photon density matrix as computed by Glauber photode- 
tection theory. This entanglement is considerable even 
for modest squeezing and highly robust against losses, 
making ~1000-partite entanglement observable with ex- 
isting methods. 

This work was inspired by pP. We thank P. Hyllus, J. 
K. Korbicz and G. Messin for helpful discussions. Sup- 
ported by the Spanish MINECO project MAGO (Ref. 
FIS201 1-23520) and by the European Research Council 
project AQUMET. 

Appendix I: Large-J S0rensen-M0lnier-Hyllus theory - 
For a given system size J, minimization of {Jh) = 
\J^ - fiJz) over spin states gives a state, parametrized 
by ^, with scaled noise V = {Jx)/J and polarization 
Z = (Jz) /J. These parametrized values describe a curve 
at the boundary between possible and impossible {V,Z) 
values for that J. k is the minimum J for which a given 
(V,Z) is possible. 

For highly polarized states, we can approx- 
imate the spin raising operator J+|J, rn) 
yJJ{J + 1) - m{m + 1) I J, m + 1) as J+ « 6\/2 J + 1 where 
h is the annihilation operator for the defect Sz = J - Jz, 
i.e., ii\Sz) = \J,J-Sz), then b\Sz) = VS^\Sz- I) ■ This 
geometrical change neglects terms of order \/V where 

v^sz/jm- 

Using Jj = ( J+ + J-)2/4, J_ = 4, we find Jh ~ (2J + 
l){b + b^)^/4 + ^b^b is quadratic in b, b^ and thus reducible 
to Jh = Ac'c + i? by a Bogoliubov transformation b -* 
ccoshr + c^ sinhr. The minimum state \(j>) is squeezed 
vacuum satisfying c|(/)) = 0, and as such has 2(jJ)/J = 

1 + 2Sz - 2^Sz{l + 5z) ^. Allowed states thus satisfy 
2V > 1 + 2JV - 2yJjV{\ + JV) or A: « J(l - 2Vfl{%V8z) 
for V < 1/2. 

With the squeezed thermal state, 5'o,S'2;,var(S'j,) play 
the role of J,Jz,vai{Jx), respectively and we have Sz = 
So - Sx = {us + nth + 2nsnt h)/2 and V = var(S'z)/5o = 
|nc(l + 2ns - 2^ns{ns + 1))(1 + 2nth)/("c + rig + nth + 
2nsnth)- In the macroscopic limit ric ->• oo, fc oc (5'o) (see 
Figure [2]). 

Appendix II: correlation functions - From [39] pp. 92-94, 
we can evaluate normally-ordered correlations in terms 



of the Wigner distribution W as 

Em,n = ((0^)"^") = fdxdpW{x,p)0^^^{x,p), (9) 



Om,7i = hm 



di3 + 



P* 



]1- 



dp* + 



o-in ifi* a+i[ja* 



,(10) 



a = X + ip. This evaluates to the efficiently computed 

r\ - ry-n-m ym 1^120^ q _ q* where 

T,„,„,a = (-2)''-™n!m!/[a!(a + n - m)!(m - a)!] . (11) 
Using the Wigner distribution of a squeezed thermal state 

2x^ , , 2p^A'^ 



2 
^2 



W{x,p) = -^^^ exp[-;^— ^] exp[ ^^^], 



rM2 



y2 



(12) 



and recalling A = \/l + n^ + \/^, T^ = 1 + 2nth , we find 

^m,7i ~ 2^a=0 7n,n,a-'-7n,n,a wncre 

Irn,n..a= f dxdpW{x,p){x + ipr-"\x^+p^r. (13) 

Defining t = {n- m)/2, these evaluate to 

r _ \^ \^ I -^^ II "^ \(~1)^_ A2(2u+2b-t-a)rp2{t+a) 

x{2{t -u + a-b)- 1)!!(2(m + b) - 1)!! (14) 

for n- m even, and zero otherwise. 

Appendix III: Reduced density matrix - We calculate 
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i2.N) 

Pl,P2),{p[,p'2\ 



E (4i«: 



t .t .t 



. . (Jirir, (Jilt 



{p3...pr,}e{H,V} 



(15) 

with respect to the state in Eq. (fTl) . Defining r, s as 
the number of V polarizations in {pi,P2} and {pi,P2}, 
respectively, we find 



K 



i2,N) 

P1;P2},{P[,P'2\ 



rn=0 \ / 



N-m-r N-m-s 

"if 



xE, 



r+m^s+m- 



(16) 



Calculations in Mathematica are available as an ancil- 
lary file, "ExtremeSpinSqueezingForPhotons.nb" 
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